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ABSTRACT 

Ribosome biogenesis in yeast requires 75 small nucleolar RNAs (snoRNAs) and a myriad of cofactors for processing, modification, 
and folding of the ribosomal RNAs (rRNAs). For the 19 RNA helicases implicated in ribosome synthesis, their sites of action and 
molecular functions have largely remained unknown. Here, we have used UV cross-linking and analysis of cDNA (CRAC) to reveal 
the pre-rRNA binding sites of the RNA helicase Rok1, which is involved in early small subunit biogenesis. Several contact sites 
were identified in the 18S rRNA sequence, which interestingly all cluster in the "foot" region of the small ribosomal subunit. 
These include a major binding site in the eukaryotic expansion segment ES6, where Rokt is required for release of the snR30 
snoRNA. Rokl directly contacts snR30 and other snoRNAs required for pre-rRNA processing. Using cross-linking, ligation and 
sequencing of hybrids (CLASH) we identified several novel pre-rRNA base-pairing sites for the snoRNAs snR30, snRIO, U3, 
and U14, which cluster in the expansion segments of the 18S rRNA. Our data suggest that these snoRNAs bridge interactions 
between the expansion segments, thereby forming an extensive interaction network that likely promotes pre-rRNA maturation 
and folding in early pre-ribosomal complexes and establishes long-range rRNA interactions during ribosome synthesis. 
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INTRODUCTION 

The synthesis of cytoplasmic ribosomes in Saccharomyces cer- 
evisiae is initiated by RNA polymerase I-mediated transcrip- 
tion of the 35S ribosomal RNA precursor (pre-rRNA), 
which contains the sequences of the mature 18S, 5.8S, and 
25S rRNAs (Fig. lA; Henras et al. 2008; Thomson et al. 
2013; Woolford and Baserga 2013). This coincides with the 
recruitment of early ribosomal proteins to the nascent tran- 
script and formation of the small subunit (SSU) processome 
from pre-assembled subcomplexes and individual biogenesis 
cofactors (for review, see Phipps et al. 2011). Initial cleavages 
(at the sites Aq, Ai , A2) (Fig. 1 A) in a complex sequence of pre- 
rRNA processing and modification events lead to the separa- 
tion of the biogenesis pathways of the large (LSU, 60S) and 
small (SSU, 40S) ribosomal subunits. In addition to 75 small 
nucleolar RNA-protein complexes (snoRNPs), which medi- 
ate both early pre-rRNA cleavages as well as most modi- 
fication events, >200 protein cofactors are involved in 
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ribosome production, among them 19 RNA helicases (Wat- 
kins and Bohnsack 2012; Martin et al. 2013; Rodriguez- 
Galan et al. 2013; Thomson et al. 2013). These RNA heh- 
cases have been proposed to act either in the structural remod- 
eling of pre-ribosomal intermediates or in the unwinding 
of snoRNA-pre-rRNA base-pairing (Ripmaster et al. 1992; 
Venema and ToUervey 1999; Martin et al. 2013). Indeed, sev- 
eral RNA helicases are required for the release of individu- 
al snoRNAs from pre-ribosomes (Kos and ToUervey 2005; 
Liang and Fournier 2006; Bohnsack et al. 2008), and depletion 
of the RNA helicase Prp43 leads to the accumulation of several 
snoRNAs on pre-60S complexes (Bohnsack et al. 2009). How- 
ever, the binding sites and molecular functions of most pre- 
ribosomal RNA helicases have remained elusive so far. 

Based on conserved sequence elements (boxes), snoRNPs 
have been grouped into the box C/D snoRNPs that mostly 
catalyze 2'-0-methylation and the box H/ACA snoRNPs 
that can mediate pseudouridinylation of the rRNA (Watkins 
and Bohnsack 2012). Besides snoRNP functions in rRNA 

©2014 Martin et al. This article, published in RNA, is available under a 
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FIGURE 1. (Legend on next page) 
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modification, several snoRNAs (U3/snR17, snR30, U14/ 
snR128, and snRlO) form base-pairing interactions with 
the pre-rRNA that are essential for the early events in pre- 
rRNA processing (see, for example, Morrissey and ToUervey 
1993; Enright et al. 1996; Mishra and Eliceiri 1997). The 
snoRNA-pre-rRNA base-pairing directs the enzyme of the 
snoRNP complex to the specific target nucleotide and the in- 
teracting regions between the RNAs were recently found to be 
more extensive than previously thought (van Nues et al. 
2011). Interestingly, modified residues cluster in the catalytic 
regions of the ribosome, such as the peptidyltransferase cen- 
ter, and neighboring modifications often require overlapping 
base-pairing sites of the snoRNAs that guide these modifica- 
tions. This could imply that snoRNAs are recruited and func- 
tion sequentially but it also raises the possibility that several 
snoRNAs may compete for base-pairing to common rRNA 
residues. 

Here, we identify the binding sites of the DEAD-box RNA 
helicase Rokl on pre-ribosomal RNA by UV cross-linking 
and analysis of cDNA (CRAC). Consistent with its role in 
small subunit biogenesis, we show that Rokl interacts with 
several distinct sites in the 18S rRNA, which cluster at the 
"foot" of the small ribosomal subunit. The major Rokl bind- 
ing site is located in the eukaryotic expansion segment ES6, 
where Rokl is required for the release of the essential 
snoRNA snR30. Rokl also strongly associates with several 
snoRNAs, and by analysis of sequence hybrids we reveal a 
number of novel pre-rRNA base-pairing sites for the 
snoRNAs U3, U14, and snRlO on the 18S rRNA. Together, 
our data imply that, possibly along with various scaffold pro- 
teins, an extensive network of long-range snoRNA-pre-rRNA 
interactions might mediate the folding and structural remod- 
eling of pre-ribosomal complexes. 

RESULTS AND DISCUSSION 

Rokl cross-links to specific sites in 18S ribosomal RNA 

To identify RNA interaction sites of the yeast RNA helicase 
Rokl we made use of the CRAC method (Granneman et al. 
2009; Bohnsack et al. 2012). Here, plasmid-encoded HTP- 
tagged Rokl (N-HTP-Rokl or Rokl-C-HTP) was expressed 
and confirmed to rescue the growth defects caused by deple- 



tion of endogenous Rokl (Bohnsack et al. 2008), demon- 
strating that tagged Rokl is able to functionally replace the 
genomically encoded Rokl. UV radiation was applied to cells 
pelleted from cultures (in vivo cross-linking; Bohnsack et al. 
2009) or to actively growing cultures in medium (in culture 
cross-linking; Bohnsack et al. 2012). The in culturo cross- 
linking method was established to reduce the manipulation 
of cells (e.g., by centrifugation, resuspension at high cell den- 
sity) before cross-linking, in order to obtain cross-links un- 
der as near steady state conditions as possible. Rokl-RNA 
complexes were then isolated by tandem-affinity purification 
where the second step was performed under denaturing con- 
ditions to ensure that RNAs covalently cross-linked to, and 
therefore directly contacting, Rokl would be isolated. 
Sequence libraries of bound RNA fragments were generated 
as previously described and analyzed by lUumina deep se- 
quencing (Bohnsack et al. 2012). The sequence reads ob- 
tained were mapped on the yeast genome as previously 
reported (Bohnsack et al. 2009). 

The distribution of Rokl cross-linking sites on the primary 
35S pre-rRNA transcript is shown in Figure 1B,D (in vivo 
cross-linking) and Figure IC (in culturo cross-linking). 
The highest hit density (peaks) was found in the 18S rRNA 
sequence, consistent with the presence of Rokl in early 90S 
and pre-40S pre-ribosomes and its role in the initial pre- 
rRNA processing events required for SSU formation (Aq, 
Al, and Aj) (Fig. lA; Jakob et al. 2012; Lebaron et al. 2013; 
Segerstolpe et al. 2013). All three approaches showed a very 
similar distribution of hits across the pre-rRNA transcript, 
confirming that Rokl cross-linking to these sites is unaffected 
by the position of the tag or the cross-linking method used. 
Within the 18S rRNA sequence, several prominent peaks be- 
tween positions 800 and 1600 were readily detectable (zooms 
in Fig. IB-D). Point mutations, introduced by the reverse 
transcriptase at nucleotides covalently linked to amino acid 
residues during UV cross-linking, were detected within all 
putative Rokl binding sites in the 18S rRNA sequence (Fig. 
IB-D zooms/lower graph), indicating direct RNA-protein 
contacts at these positions. Rokl is the second RNA helicase, 
besides Prp43 (Bohnsack et al. 2009), for which binding sites 
were identified on yeast ribosomal RNA, and these two pro- 
teins show very different cross-linking patterns, indicating 
the specificity of the sites identified. 



FIGURE 1 . The Rokl cross-linking sites on 18S rRNA cluster in the 3D structure of the small ribosomal subunit. (A) Schematic overview of the key 
steps in pre-rRNA processing pathway in S. cerevisiae. ETS, external transcribed spacer; ITS, internal transcribed spacer; Nu, nucleus; Cy, cytoplasm. 
(B-D) Cells expressing plasmid encoded N-terminally {B,C) or C-terminaUy (D) HTP-tagged Rokl were depleted of endogenous Rokl and UV cross- 
linked in vivo {B,D) or in culturo (C). Rokl -containing complexes were purified, cross-linked RNAs were trimmed, and linkers ligated, followed by 
RT-PCR and Illumina deep sequencing. Obtained sequence reads were mapped on the yeast genome and the results for RDN37 encoding the 35S pre- 
rRNA are plotted as number of hits per million reads for each nucleotide. Relative positions of the mature 18S, 5.8S, and 25S rRNA sequences are 
indicated below. A magnification of the major Rokl cross-linking sites is given. Mutations arising during reverse transcription at cross-linking sites 
are shown in the lower panel of the zoom. The asterisk indicates a peak that is also present in untagged control samples. (£) Rokl cross-linking sites (N- 
HTP in culturo) with hit values above a threshold of 5% of the highest peak on 35S rRNA are visualized as colored circles on the secondary structure of 
the mature 18S rRNA. Colors indicate peak height with the highest peak (100%) shown in red and lesser peaks shown in yellow. Nucleotides marked in 
green indicate a peak that was just above/below the 5% threshold in the different data sets. (F) Rokl cross-linking sites (N-HTP in culturo) mapped on 
the 3D structure of the small ribosomal subunit. Peaks on rRNA are colored as in E and the location of cross-linking sites is indicated. H; helix; ES, 
expansion segment. Ribosomal proteins Rpsl3 and Rpsl4 are shown in purple. 
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Mapping the Rokl cross-linking sites onto the 2D structure 
of the mature 18S rRNA revealed that the main cross-linking 
sites of Rokl are found in the eukaryotic expansion segment 6 
(ES6) on helices (H) 1 and 3 (ES6H1 andES6H3, respectively) 
(Fig. IE). The other cross-linking sites of Rokl are located in 
Hll, eukaryotic expansion segment 3 (ES3), H9 that im- 
mediately precedes ES3, and H27 of 18S rRNA (Fig. IE). In 
addition, cross-linking between Rokl and sequences in eu- 
karyotic expansion segment 7 (ES7) was observed in all exper- 
iments with peaks just above or just below the 5% threshold. 
Although the Rokl cross-linking sites are distributed across 
the 5' end of the 1 8S rRNA sequence, the formation of second- 
ary and tertiary structures during ribosome synthesis can 
bring such distant regions of the pre-rRNA transcript togeth- 
er. We therefore mapped the cross-linking sites of Rokl onto 
the 3D structure of the yeast 40S ribosomal subunit (Fig. IF; 
Ben-Shem et al. 2011). Remarkably, all of the Rokl cross-link- 
ing sites besides that in H27, which is only just above the 
threshold, clustered on one face of the "foot" region of the 
small ribosomal subunit (SSU). Furthermore, two ribosomal 
proteins, Rpsl4 and Rpsl3, have previously been shown to be 
required for Rokl association with pre-ribosomes (Jakob et al. 
2012). Interestingly, these proteins are localized adjacent to 
newly identified Rokl cross-linking sites (Fig. IF). 

Taken together, the identification of Rokl-pre-rRNA 
cross-linking sites in the IBS rRNA sequence by our CRAC 
analyses is consistent with the role of this helicase in SSU 
maturation. Furthermore, the close proximity of the distinct 
Rokl cross-linking sites in the 3D structure of the SSU might 
indicate that they are derived from a single Rokl binding site 
on the pre-ribosome. 
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FIGURE 2. DMS structure probing confirms the interaction site of 
Rokl in ES6H3 of the 18S rRNA sequence. (A) Extracts from cells ex- 
pressing HA-tagged Rokl under tetracycline control and Enpl-HTP 
were incubated with IgG sepharose and the eluate was analyzed by 
SDS-PAGE and Western blotting using antibodies against the HA-tag 
(Rokl), HTP-tag (Enpl), or endogenous GAPDH. (B) RNA from the 
Enpl-HTP pulldown experiments described in A was isolated and ana- 
lyzed by Northern blotting using a probe hybridizing to the 5' end of the 
internal transcribed spacer 1 (ITSl). Pre-rRNAs (indicated to the left) 
were visualized using a phosphorimager, and mature 18S and 25S 
rRNAs were detected by methylene blue staining. (C) Cells expressing 
Rokl or depleted of Rokl were treated with DMS in vivo and RNA 
that co-precipitated with Enpl-HTP (as in A) was isolated. RNA was an- 
alyzed by primer extension using a primer annealing downstream from 
ES6H3 to detect DMS-modified nucleotides (numbers indicate residues 
within the 18S rRNA sequence). Labelled RNA fragments were separated 
by denaturing polyacrylamide gel electrophoresis and visualized by 
phosphorimaging. The position of the Rokl cross-linking site is indicat- 
ed on the left. 



Structure probing confirms Rokl binding to ES6H3 

To validate the pre-rRNA cross-linking sites of Rokl identi- 
fied here by CRAC we performed in vivo chemical probing 
(Wells et al. 2000), which allows us to detect structural chang- 
es within the predicted pre-rRNA interaction site upon Rokl 
binding. Treatment of actively growing cells with DMS leads 
to methylation of accessible nucleotide residues (predomi- 
nantly adenosines) and such modifications can be detected 
as primer extension stops. As the Rokl cross-linking site is lo- 
cated within the mature 18S rRNA sequence, we isolated pre- 
rRNAs from Rokl -containing pre-ribosomal complexes to 
enable detection of changes in DMS-induced modification 
caused by Rokl binding. Enpl is an SSU biogenesis factor 
that is associated with both early and late pre-ribosomal com- 
plexes making it an ideal factor through which to isolate 
Rokl -containing particles. A yeast strain expressing genomi- 
cally tagged Enpl-C-HTP was therefore used in pulldown ex- 
periments, and the presence of Rokl as well as the enrichment 
of the 35S and 20S pre-rRNAs and lack of mature 18S rRNA 
were confirmed by Western and Northern blotting, respec- 
tively (Fig. 2A,B). RNA associated with Enpl particles was 
isolated from DMS-treated cells either expressing, or depleted 



of, Rokl. Reverse transcription using a primer located down- 
stream from the Rokl cross-linking site in ES6Fi3 showed 
that the DMS modification of multiple nucleotides (such as 
A803, A809, A811, and A850) was unaffected by depletion 
of Rokl (Fig. 2C). Excitingly, the accessibility of residues lo- 
cated within the Rokl binding site (A817, G823, A829, and 
U832) was significantly increased (Fig. 2C). These data con- 
firm that the Rokl cross-linking site identified by CRAC rep- 
resents bona fide interactions with the 18S rRNA sequence. 

Interestingly, ES6H3 and sequences nearby were also re- 
cently identified as cross-linking sites of two other cofactors, 
Rrp5 and Rrp7, which are also involved in the early steps of 
ribosome biogenesis (Lebaron et al. 2013; Lin et al. 2013). 
While Rrp7 is recruited to the 90S complex as part of a pre-as- 
sembled module, the UTP-C complex, Rrp5 has further been 
proposed to function as a central scaffold, coordinating pre- 
ribosome assembly. Consistent with their adjacent pre- 
rRNA binding sites, Rrp5 has been shown to directly interact 
with Rokl. Furthermore, Rrp5 is suggested to modulate the 
activity of Rokl and to confer target specificity on the RNA 
helicase (Garcia et al. 2012; Young et al. 2013). In the case 
of Rrp7, cross-linking sites in ES6Fi3 and ES7 that are 
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overlapping with those suggested here by our Rokl CRAC 
analysis were previously identified (Lin et al. 20 1 3 ) . Taken to- 
gether, these data could suggest that Rokl and Rrp7 may asso- 
ciate sequentially with pre-ribosomal complexes. 

Rok1 directly interacts with snR30 

We have previously shown that Rokl is required for the re- 
lease of the snR30 snoRNA from pre-ribosomal complexes 
(Bohnsack et al. 2008). Since snR30 base-pairs with the 18S 
rRNA sequence in ES6 (at sites rml and rm2), this pro- 
vides independent evidence for Rokl binding and function- 
ing in this region of the pre-ribosome. Furthermore, the 
overlapping binding sites of Rokl and snR30 may indicate a 
direct role for Rokl in regulating the association of snR30. 
Therefore, to gain more insight into the interactions of 
Rokl with this snoRNA-pre-rRNA duplex, we first ascer- 
tained, by analysis of our CRAC data, that snR30 strongly 
cross-links to Rokl. We then analyzed the distribution of 
the Rokl cross-linking sites along this snoRNA and found 
several (and corresponding point mutations in the sequence 
reads) along the length of the snoRNA (Fig. 3A) . Interestingly, 
the majority of peaks were localized at the 3' end of the 
snoRNA where the functionally essential elements of this 
snoRNA are found (Atzorn et al. 2004; Fayet-Lebaron et al. 
2009). Significant cross-linking of Rokl to the internal hair- 
pin that is structurally important for snR30 function was ob- 
served. The m2 motif of snR30 was also identified as a 
prominent Rokl interaction site. This correlates well with 
our identification of the corresponding IBS rRNA sequence, 
rm2, as a Rokl binding site (Fig. 3B) and could suggest that 
Rokl mediates the release of snR30 by local strand unwind- 
ing. When bound to snR30, ES6H3 is suggested to form a 
hairpin structure (Fig. 3B; Fayet-Lebaron et al. 2009). Since 
Rokl is proposed to unwind RNA duplexes with a single 
stranded extension (Garcia et al. 2012), it is possible that 
Rokl mediates unwinding of the snR30-ES6H3 duplex 
from a binding site in the loop of this hairpin. Furthermore, 
the cross-linking sites of Rokl and Rrp5 on snR30 are located 
in close proximity, which is in line with the adjacent binding 
sites of these two proteins in ES6H3 and ES6H2 of IBS rRNA, 
respectively. 

Due to its base-pairing interactions with ES6H3, the re- 
lease of snR30 from pre-ribosomal particles is a prerequisite 
for formation of the interactions between ES6H3 and ES3 
that are found in the mature ribosome (Alkemar and 
Nygard 2003; Ben-Shem et al. 2011). As the Rokl helicase 
contacts all three RNA sequences involved and is required 
for pre-ribosomal release of snR30, it may function as a key 
regulator of this structural rearrangement. Furthermore, 
both ATP-dependent unwinding and ADP-dependent an- 
nealing functions of Rokl have been described (Garcia 
et al. 2012; Young et al. 2013), possibly implicating this heli- 
case in formation of base-pairing interactions between ES6 
and ES3 after the release of snR30 from ES6. 



Novel snR30 base-pairing sites on the 1 8S rRNA 

It has recently become clear that the extent of interactions be- 
tween snoRNAs and rRNA has been significantly underesti- 
mated; additional base-pairing sequences of snoRNAs 
involved in both rRNA modification and processing have 
been identified (Kudla et al. 2011; van Nues et al. 2011). We 
therefore speculated that snR30 might also form additional 
base-pairing interactions with the IBS rRNA. To investigate 
this hypothesis, we employed two complementary approach- 
es. Firstly, computational prediction of snR30 base-pairing 
sites in the IBS rRNA sequence was performed. Secondly, us- 
ing the cross-linking, ligation and sequencing of hybrids 
(CLASH) (Kudla et al. 2011) method we analyzed chimeric 
reads, which arise due to ligation, and thus hetero-concate- 
merization, of cross-linked RNA fragments or by strand 
switching of the reverse transcriptase on cross-linked RNA 
duplexes. Base-pairing prediction has the advantage of pro- 
viding an unbiased and sensitive search, while chimeric se- 
quences in CRAC data provide experimental evidence of 
snoRNA-pre-rRNA interactions at the binding sites of pro- 
teins of interest. 

For the prediction of snR30-lBS rRNA hybrids, the se- 
quences of both RNAs were divided into short, overlapping 
fragments and screened for potential base-pairing interac- 
tions using scripts including the RNAduplex algorithm 
(ViennaRNA package) (Lorenz et al. 2011), and only hits 
above a minimum stability threshold were considered signif- 
icant (see Materials and Methods). This approach identified 
three potential base-pairing sites between snR30 and IBS 
rRNA while, excitingly, CLASH analysis of the deep sequenc- 
ing data derived from our Rokl CRAC experiments revealed 
hybrids corresponding not only to the known snR30 site in 
ES6H3 but also additional putative interactions sites in 
ES6H1, ES7, H15, H22, and H44 of IBS rRNA (Fig. 3C,D). 
The sites identified by these two independent approaches 
were overlaid graphically according to their positions on 
the linear RNA sequences of snR30 relative to IBS (Fig. 3C) 
and also mapped onto a secondary structure map of IBS 
rRNA (Fig. 3D). In addition to the known snR30 base-pair- 
ing sites in ES6H3 (CI), this excitingly revealed two putative 
new snR30 interaction sites (C2 and C3), supported by the 
identification of overlapping CLASH hybrids and predicted 
base-pairing sites on both IBS rRNA and snR30 (Fig. 3C, 
D). Firstly, we find a prediction-CLASH hybrid pair that 
overlap both on snR30 and in ES6H1 of the IBS rRNA 
(C2). Secondly, our base-pairing prediction analysis also 
highlighted a very stable, 19-nt sequence in ES7 as a potential 
snR30 base-pairing site that is supported by the identification 
of an overlapping CLASH hybrid (C3). This potential inter- 
action is of particular interest because it involves both the 
snR30 sequence most frequently cross-linked by Rokl and 
also a putative Rokl contact site in the IBS rRNA sequence. 
Furthermore, this potential base-pairing site overlaps with 
a recently identified interaction site of the U3 snoRNA in 
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FIGURE 3. snR30 directly interacts with Rokl and has additional base-pairing sites in the 18S rRNA sequence. (A) Rokl-N-HTP reads that mapped 
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ES7 (Kudla et al. 2011). In addition, base-pairing prediction 
analysis of Schizosaccharomyces pombe snR30 and 18S rRNA 
also identified a likely interaction (minimum free energy 
— 16 kcal/mol) between ES7 of the 18S rRNA sequence (nu- 



cleotides 1051-1070) and a region of the snR30 internal hair- 
pin (nucleotides 226-245) also identified in S. cerevisiae. This 
implies that this potential interaction site of snR30 is likely 
conserved in other fungi. 
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Rok1 -associated snoRNAs form multiple interactions 
with the 1 8S rRNA 

Besides snR30, pre-rRNA interaction sites of other essential 
snoRNAs have also been described; U3 interacts with 18S 
rRNA sequences that ultimately form the central pseudoknot 
as well as with sites in the 5' external transcribed spacer (ETS), 
while U14 base-pairs at its modification site in H13 and also 
in H6a. Less is known about the base-pairing sites of snRlO 
as the short stretches of base-pairing interactions that box 
H/ACA snoRNAs form with the pre-rRNA often make their 
bioinformatic identification challenging. Although a binding 
site in the 5' ETS has been reported (Liang et al. 2010), poten- 
tial interactions of snRlO with the 18S rRNA remain elusive. 
Interestingly, a significant enrichment of U14, U3a, and to a 
lesser extent, snRlO sequences was found in our Rokl 
CRAC analysis, implying that, as well as snR30, Rokl might 
interact with these snoRNAs (Fig. 4A). This raised the possi- 
bility that, like snR30, these snoRNAs may base-pair with the 
pre-rRNA at Rokl binding sites. 

Although Rokl did not appear to interact with snRlO as 
strongly as with snR30, U3a, or U14, genetic interactions be- 
tween snRlO and Rokl have previously been reported indicat- 
ing a functional link (Venema et al. 1997). Analysis of the 
distribution of Rokl hits along snRlO revealed two peaks: 
one adjacent to an essential 7-nt sequence (Liang et al. 
2010) and the Rrp5-binding site (Lebaron et al. 2013) and a 
second at the 3' end of the snoRNA where other functionally 
important motifs such as the ACA box are found (Fig. 4B). 
This indicates that Rokl may be a regulator of snRlO func- 
tion, so we used CLASH to determine if snRlO base-pairs 
with 18S rRNA at any of the Rokl binding sites. This uncov- 
ered two potential snRlO interaction sites in ES6H1 and 
ES6H3 (Fig. 4C); however, CLASH hybrids between ES6H3 
and snRlO were not recovered as frequently as those arising 
between ES6H3 and snR30 described above. Interestingly, 
snRlO was recently reported to interact strongly with Rrp7 
and, based on the binding of Rrp7 to ES7, an snRlO base- 
pairing site in this region was proposed (Lin et al. 2013). 
Since Rrp7, like Rokl, also binds to ES6H3, it is tempting 
to speculate that snRlO, possibly regulated by Rokl or 
Rrp7, bridges interactions between ES6H3 and ES7 in pre-ri- 
bosomal complexes. 

For U14, the CLASH analysis identified multiple hybrids 
containing the previously reported U14 base-pairing and 
modification site in H13 (Fig. 4D). However, additional hy- 
brid sequences overlapping with the Rokl cross-linking sites 
in H9 and Hll were also found. Rrp5 interacts with the IBS 
rRNA sequence in between these two binding sites (Lebaron 
et al. 2013), again suggesting an interplay of these factors 
within these regions of the pre-ribosome. 

We also detected hybrids between the U3 snoRNA and 
ES6H3, ES6H1, ES3, and H9 of the 18S rRNA sequence 
(Fig. 4E). Depletion of Rokl has been shown to cause some 
accumulation of U3 in pre-ribosomal complexes (Bohnsack 



et al. 2008), which might imply a functional relevance of 
these interactions. It is important to note that the known 
U3-binding sites in the 18S rRNA sequence or the 5' ETS 
do not overlap with the major Rokl binding sites. We were 
therefore unlikely to identify these sequences in our analysis. 
Interestingly, the discovery of potential U3 base-pairing sites 
in both ES3 and ES6H3/1 raises the possibility that the U3 
snoRNA bridges long-range interactions between these se- 
quences in pre-ribosomes. These interactions, possibly regu- 
lated by Rokl, might contribute to folding of the pre-rRNA 
and help establish the base-pairing between ES3 and ES6, 
which is found in mature ribosomes. 

RNA folding and compaction of early pre-ribosomal 
complexes into "terminal knobs" can be observed in Miller 
spreads (Osheim et al. 2004), and recently Rrp5 was proposed 
to function as a scaffold, coordinating pre-rRNA cleavages 
and assembly of pre-ribosomes (Lebaron et al. 2013). In con- 
cert with this proteinaceous framework, snoRNAs are emerg- 
ing as additional players in establishing structures within 
pre-ribosomal complexes that are essential for maturation. 
On the one hand, interactions between snoRNAs and pre- 
rRNA are likely important for maintaining an open structure 
of early pre-ribosomes that is required for rRNA modifica- 
tion by other snoRNPs and the association/action of many 
biogenesis factors. On the other hand, snoRNAs have im- 
portant functions in coordinating pre-rRNA cleavages, likely 
acting in local RNA folding and bringing together distant re- 
gions of the transcript. More specifically, snR30, snRlO, and 
U3a base-pairing sites appear to cluster in both ES7 and ES6 
(Kudla et al. 2011; Lin et al. 2013), suggesting that these re- 
gions might be tethered together in pre-ribosomal complex- 
es. Interestingly, some of the new base-pairing sites of these 
snoRNAs are adjacent or even overlapping, perhaps indicat- 
ing sequential association of the snoRNAs. Indeed, it is pos- 
sible to imagine a scenario in which ES6 is first structurally 
linked to ES7 by several snoRNAs and that subsequent 
snoRNA-mediated action brings together ES6 and ES3 to 
establish the structure found in mature ribosomes. We there- 
fore propose a snoRNA-pre-rRNA-protein interaction net- 
work that orchestrates the remodeling of pre-ribosomal 
complexes (Fig. 4F). As an active ATP-dependent RNA heli- 
case Rokl is likely to provide the catalytic activity required to 
drive some of the early transitions in the structure of the pre- 
ribosome. Interestingly, Rrp5 is one of the few ribosome bio- 
genesis factors that is implicated in the maturation and in fact 
contacts rRNAs of both ribosomal subunits (Venema and 
ToUervey 1996; Lebaron et al. 2013). This raises the possibil- 
ity that the snoRNAs involved in the pre-SSU scaffold might 
not only coordinate rearrangement of the early pre-SSU but 
also form contacts to pre-rRNAs of the LSU, thereby cou- 
pling maturation of both ribosomal subunits. 

In prokaryotes, ribosomal proteins themselves largely 
coordinate maturation and structuring of the ribosomal 
subunits while the greater size, complexity, and number of 
rRNA modifications of eukaryotic ribosomes appear to 
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FIGURE 4. CLASH hybrid analysis identifies novel 18S rRNA base-pairing sites for snRlO, U14, and U3. (A) Heatmap showing the relative cross- 
linking of Rokl to snoRNAs of the SSU. (B) Rokl-N-HTP reads that mapped to snRlO are shown as the number of hits for each nucleotide per million 
mapped reads. The asterisk indicates a peak also found in control samples. Mutations arising during reverse transcription at cross-linked nucleotides 
are shown in the lower panel. Black bars below indicate positions of conserved sequence elements and the hash indicates the position of the functionally 
important 7-nt sequence of snRlO. The bracket indicates the position of the previously identified Rrp5 interaction site on snRlO. {C-E) The positions 
of identified CLASH hybrids on the secondary structure of 18S rRNA are shown for snRlO (C), U14 (D), and U3a (E). (F) Interaction network of 
relevant 18S rRNA regions (left), snoRNAs (right), and proteins (Rokl, Rrp5, and Rrp7; center). Dotted line indicates direct protein-protein inter- 
action, solid lines mark RNA-protein interactions identified by CRAC, and dashed lines denote snoRNA-pre-rRNA interactions detected by 
CLASH. The asterisk indicates that Rrp5 interacts with an 18S rRNA sequence adjacent to, but not within, ES6H3. 



necessitate the evolution of the large number of ribosome 
biogenesis cofactors and the formation of a scaffold of 
snoRNAs and proteins to regulate ribosome assembly. One 
of the major differences between prokaryotic and eukaryotic 
rRNAs is the presence of eukaryotic expansion segments, 
and intriguingly most of the newly identified snoRNA- 
pre-rRNA interactions revealed here cluster in these re- 
gions, suggesting that they may act as important anchor 
points from which eukaryotic pre-ribosome remodeling is 
coordinated. 



MATERIALS AND METHODS 

UV cross-linking and analysis of cDNA (CRAC) 

For CRAC experiments plasmids based on pRS415 that enable ex- 
pression of either N- or C-terminaUy Hise-Tev-ProteinA (HTP)- 
tagged Rokl were transfected into a Rokl depletion strain, 
YMB 146 (Bohnsack et al. 2008). Where necessary, Rokl was depleted 
by addition of 20 |-ig/mL doxcycline for 9 h. The growth rates of 
all strains were determined as previously described (Bohnsack et 
al. 2008). CRAC experiments and analysis of the deep sequencing 
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was performed as previously described (Bohnsack et al. 2009; 
Granneman et al. 2009; Bohnsack et al. 2012). In brief, UV cross- 
linking was performed in living cells either in culture medium (in 
culturo) or pelleted and resuspended in small volume (in vivo), 
and Rokl -containing complexes were isolated first on IgG sepharose, 
eluted with TEV protease, and purified on Nickel-NTA. After RNase 
trimming and linker ligation, fragments were amplified by RT-PCR 
and sent for Ulumina deep sequencing. Mapping of cross-linking 
data onto the 18S rRNA secondary structure (Petrov et al. 2013) 
was performed by using Python (Version 2.7) scripts and vector- 
graphics from http://apollo.chemistry.gatech.edu/RibosomeGallery. 
Mapping of cross-linking data on the 3D structure of the SSU 
(Ben-Shem et al. 2011) was performed using Python scripts and 
Pymol. The height of all peaks relative to the highest peak was calcu- 
lated and they are mapped in a color gradient of yellow (5%) to red 
(100%). 

Isolation of Enpl -containing pre-ribosomal 
complexes 

For purification of pre-ribosomal complexes chromosomally en- 
coded Enpl was C-terminally HTP-tagged in the YMB146 strain. 
Exponentially growing cells were harvested and lysed by grinding 
in liquid nitrogen. Enpl-HTP-containing pre-ribosomal complexes 
were isolated in a one-step purification on IgG sepharose before iso- 
lation of co-precipitated RNA using phenol-chloroform. Northern 
blot analysis was performed using a probe hybridizing to the 
5' end of ITSl (5'-CGGTTTTAATTGTCCTA-3'). For analysis of 
the protein composition of purified pre-ribosomal particles, a sam- 
ple of the pulldown eluate was separated by SDS-PAGE and analyzed 
by Western blotting using anti-HA (Rokl; Sigma- Aldrich), anti- 
PAP (Enpl; Sigma- Aldrich), and anti-GAPDH (Abeam) antibodies. 

In vivo chemical probing 

DMS structure probing using the Rokl depletion/Enpl-C-HTP 
strain was performed as previously described (Swiatkowska et al. 
2012). Briefly, exponentially growing cells were depleted of Rokl 
for 9 h before pelleting and resuspension in 25 mL YPD medium. 
DMS was added to a final concentration of 1% for 4 min at 30°C be- 
fore quenching with 0.7 M 6-Mercaptoethanol in isoamylalcohol 
and washing with water. RNA was isolated from pre-ribosomal 
complexes purified via Enpl-HTP as described above. Primer exten- 
sion using a radiolabeUed primer (5'-CGTCCTTGGCAAATGC-3') 
was performed using Superscript III reverse transcriptase (Life 
Technologies) according to the manufacturer's instructions. RNA 
was separated on a 10% denaturing, polyacrylamide sequencing 
gel and visualized using a phosphorimager. A sequencing ladder 
was prepared using the above oligonucleotide, a pAV162 plasmid 
containing the 35S rDNA as a template and the Sequenase 
Version 2.0 DNA Sequencing Kit (Affymetrix) according to the 
manufacturer's instructions. 

CLASH analysis and hybrid predictions 

CLASH analysis of hybrids formed during the ligation step of the 
CRAG procedure was performed as previously described (Kudla 
et al. 201 1; Travis et al. 2013). In brief, sequences were pre-processed 
using Flexbar (trim = 0, fdter = 0, minimum = 4, length = 17) and 



reads were mapped using BLASTALL. The computational 
snoRNA-rRNA interaction prediction algorithm was based on the 
RNAduplex software (ViennaRNA package; Lorenz et al. 2011) 
and was adapted into a Python script designed to test the base-pair- 
ing properties for 20-nt RNA fragments with a 15-nt overlap. A cut- 
off of —20 kcal/mol minimum free energy was used as a threshold to 
define stable hybrids. Clustering of experimentally obtained hybrids 
and predicted base-pairing sites was performed with Python scripts 
and matplotlib (Hunter 2007). 
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